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We theoretically study bipartite entanglement purification with neutral atoms via cavity-assistant 
interaction and linear optical elements. We focus on entanglement distillation and the recurrence 
protocol, whose performances under idealized and realistic conditions are discussed. The implemen- 
tation of these purification protocols has been tested with numerical simulations. We analyze the 
performance and stability of all required operations and emphasize that all techniques are feasible 
with current experimental technology. 

I. INTRODUCTION 

Entanglement has always been at the heart of quantum mechanics and plays a key role in quantum information 
science. Most applications of quantum information processing require maximally entangled states. Many efforts have 
been undertaken to generate entangled states in different physical systems. Recently, the controlled manipulation 
of certain systems in such a way that entangled states can be generated on demand has become possible. However, 
noise in such control operations as well as interactions with an uncontrollable environment result that the desired 
entangled states are created only with a certain and less than unit fidelity. There are several ways toprotect quantum 



information and entangled quantum states in particular including quantum error correction 



al3|,lJ,l5|,l6|,l7l 



12| and experimentally demonstrated in 
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and entanglement purification which have been derived in 
Here we will focus on the latter method and the physical implementation with neutral atoms. 

A number of different entanglement purification protocols exist, which differ in their purification range, the efficiency, 
and the number of copies of the states they operate on. In the following, we show the physical implementation of en- 
tanglement purification on neutral atoms via cavity-assistant interaction and linear optical elements. We demonstrate 
the distillation of maximally entangled states from both non-maximally and mixed entangled inputs. In Sec. H, we 
implement a positive operator valued measure (POVM) and full Bell state measurement (BSM) on atoms via cavity- 
assistant interaction and linear optical elements, which are basic tools for entanglement purification. In Sec. HI, we 
will consider bipartite entanglement distillation (also called as filtering protocols) which operates on a single copy, as 
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well as the recurrence protocol which operates on two copies simultaneously at each step in Sec. IV. With linear opti- 
cal elements, we perform a filtering process to maximize the entanglement of pure atomic entangled states. We have 
also applied our methods to initial states which are partially mixed. After filtering, the distilled states demonstrate 
certain non-local correlations, as evidenced by their violation of a form of Bell's inequality. Because the initial states 
do not have this property, they can be said to possess "hidden" non-locality. Furthermore, we propose the physical 
implementation of the recurrence protocol with the same setup. In Sec. V, the implementation of the proposals have 
been tested with numerical simulations, which included various sources of noise present in experiments, and was found 
to be robust against the influence of noise. We analyze the feasibility of the implementation of these puriflcation pro- 
tocols and found it practically based on the POVM and BSM with high fidelities in the current status of experimental 
technology. Remarkably, our protocol will enable the first-ever implementation of atomic entanglement distillation 
and purification. 

II. THEORETICAL MODEL AND BASICS TOOLS 
A. Theoretical Model 

Our source of entanglement and some basics tools which are used to implement entanglement purification include 
a POVM and BSM, which are taken as cavity-assistant interaction on neutral atoms. The atoms are trapped in a 
transverse optical lattice and thus can be taken into and out of the cavity for local operations in order to circumvent the 
requirement for individual addressing. We consider a high Q optical cavity couples photons to the certain transitions 
in atoms whose internal states are shown in Fig. [TJa. The key parameters describing a cavity QED system are the 
cavity resonance frequency Wc, the frequency tJa = ("-^c) of the atomic transition from the ground state |0) to the 
excited state |e). The other ground state |1) is decoupled due to the large hyperfine splitting. The system of the 
cavity and atoms are described by the Hamiltonian (setting h ~ 1) 

H = Ho + Hi^t + H., + H^, (2.1) 

where 

i/o - uj,a)a + ^c^a(|e)^. (e| - |0)^. (0|); 
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Hint = S3 («c (0| + 4 |0)^. (e|) , (2.2) 
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and describes the coupling of the cavity to the continuum which produces the cavity decay rate n = lUc/Q, while 
Hj describes the coupling of the atom to modes other than the cavity mode which cause the excited state to decay 
at rate 7. 

We use the single-photon pulse resonant with the bare cavity mode with the horizontal polarization h. Thus the 
input state of cavity is \h). Suppose we perform this operation in the limit with T 3> that is the pulse duration 
T is much large than the time that the photon state in the cavity. We now present a detailed theoretical model with 
one atom as an example. In the rotating wave approximation, the Hamiltonian of atom-cavity and free space is 



Hsy = - i| |e) (e| + gia^ |e) (0| + h.c.) + Aala^ 

Ludujb^ (uj) b (lj) 



(2.3) 



where A denotes the detuning of the cavity field mode Oc from the atomic transition, and b (lu) with the standard 
relation [6 (lu) , (w')] = S {lu — lu') denotes the one-dimensional free-space modes which couple to the cavity mode 
Uc- According to the quantum Langevin equation and the boundary condition of the cavity, we can deduce that 
the single-sided cavity input and output field operators (t) and 6out(i) sie connected with the cavity mode a^it) 



through the relations 
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and 



ac{t) = -i [a^it), Hsy] - (iA + ^) adt) - \/^6i„(t) 



bout{t) = hn{t) + ^/Kac{t), 



(2.4) 



(2.5) 



where the Hamiltonian Hgy shown in Eq. (2.3) describes the coherent interaction between the atom and the cavity 
mode flc. The operators bin{t) and bout{t) satisfy the following commutation relations fein(i), bj^ (t/) = S (t — t') = 
bout{t),bl^^. {t') . If the atom is in the state |1), the Hamiltonian Hsy is not active, and induces A = 0. When 
the input pulse shape changes slowly with time t compared with the cavity decay rate k, from Eqs. (2.4) and (2.5), 
we obtain bpntit) ~ !a+k/2 ^'"(^) ~ ^hn{t)- While, if the atom is in the state |0), in the case of strong coupling 
26 
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iA+K/2 

28l |. the effective detunings of two dressed cavity modes from the input pulse are A = ±g. Consider that the 



system works in the strong coupling regime g 3> k, we obtain 6out(0 ~ ^in(i)- In practice it turned out that the result 



is true even if g ~ k [17[. From the description above, we conclude that the state of the whole system of atom-cavity 
and free-space acquires the phase tt or 0, for 6out(i) ~ ~^in(i) or 6out(^) ~ &in(i), after the photon pulse has been 



reflected by the cavity. The input-output process can be characterized by 



.i.|i,/.>(i,/.|(„|0) + 5|i))|;,)^(„|0)_6|l))|/i) 



(2.6) 



where we have discarded the state of cavity since in the hmit T ^ it is always in the vacuum state 
This model has been used to implement quantum gate and entanglement generation 
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2l| . In our 



paper we will introduce for the first time how to use this model to implement a POVM on atomic state which is not 
trivial and very useful to quantum information processing such as entanglement purification. 



B. Physical Implementation of a POVM 

According to Fig. [1^, if we choose 

{a\Q)+h\l)){\h) + \v))/V2 (2.7) 



as initial state of the atom+photon system, then the photons with horizontal and vertical polarizations are "bounced" 
back from the atom-cavity system and a mirror M respectively. After leaving the cavity and passing through the first 
polarizing beam splitter (PBSl), which, as we recall, reflects vertically polarized light \v) and transmits horizontally 
polarized light the polarization of the photon is rotated by 7r/4 using half wave plate 1 (HWPl) which performs 
\h) {\h) + \v)) /^/2 and \v) {\h) - \v)) /V2, and the state 

a\Q)\h)+b\\)\v) (2.8) 

is obtained. Afterwards, the photon passes through PBS2. The polarization of the reflected photon is subsequently 
rotated by HWP2. The extra rotation allows the applied partial polarizer to operate in the basis {\h) , \v)} as 

\h)^e\h),\v)-.\v). (2.9) 

This filtered state then interferes at PBS3 with the state that has been transmitted at PBS2. The photon passes 
HWP3 and PBS4 and is then detected by single-photon detectors. If detector D2 clicks, the resulting atomic state is 
a |0) + eb\l) (not normalized). That is we implement a POVM on the atomic state as 

P=|0)(0|+e|l)(l|. (2.10) 



The partial polarizer 



22j can be realized experimentally by inserting into one path a series of slabs, tilted about 



the vertical axis by 9 (see Appendix for details). Owing to the well known polarization dependent reflectivity, the 



5 

transmitted photons are preferentially vertically polarized. In the ideal case, the vertical polarized photons are 
perfectly transmitted, while the horizontal polarized ones are partially reflected. The transmission probability for 
horizontal polarization is Ti^ ~ e, while for vertical polarization it is only Tv = 1. 

In addition, if one adds a quarter wave plate (QWP) between partial polarizer and PBS3, another POVM 

P' = |0)(0|+i£|l)(l| (2.11) 

can be easily implemented on atomic states. 



C. Physical Implementation of full BSM 

For two atoms which are trapped in different optical lattices and both interacting with the cavity, we can implement 



2l| . Here we give more detailed description 



a BSM on two atoms, which has been mentioned briefly in our recent work 
of BSM on atoms. 

First the single photon pulse enters the cavity with only atom 1 trapped in optical lattice A. After the interaction 
between atom and cavity mode, a gate operation e"^|i^'')(i.'»l jg applied on the atom and the photon pulse. Atom 2 
trapped in optical lattice B now is moved into the cavity while atom 1 is outside, and the pulse is reflected successively 
to enter the cavity again, so that the same operation is applied on atom 2 and the pulse. That is 



Yl e-|i).<i|«l'^)<'^l(ai \0), + b, \l),)ia2 \0), + b, \l),)i\h) + \v))/V2 

i=l,2 

= (ai |0>i - 6i |l)i)(a2 |0>2 - &2 11)2) \h) + (ai \0}, + b, \l),){a2 \0}, + &2 |1>2) \v) ■ (2.12) 

Finally, the photon passes through a HWP and is detected by either Dl or D2, we obtain the projection if the detectors 
Dl and D2 click respectively 

Pi = |00) (00| + 111) (111 ; P2 = 1 - Pi. (2.13) 

Performing the measurement {Pi,P2} on atoms allows one to distinguish the subspace spanned by {|$^) , I'l' )} 
and {|^^) , 1^ )}■ The measurement outcomes Pi and P2 correspond to 

Al*+>.l*->} = 1*^) + I*") (2-14) 

and 

^{|*+>,l*-)} = 1*^) + I*") ' (2-15) 
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respectively, where 1$=*=) = (|00) ± |ll))/\/2 and j^*^) = (|01) ± |10))/%/2. More generally, one can obtain non- 
destructive projections onto subspaces spanned by two arbitrary Bell states using additional single qubit unitary 
operations which allow one to permute Bell states. For instance, the application H ® H (here H notes Hadamard 
gate) consequently before and after the measurement P||$+^ j^-^i, corresponds to ^{|*+),|*+)} because of 

H®H\<^+) = \<^+) ,H ® H\^-) = (2.16) 

Obviously, using these non-destructive projections, we can achieve a full BSM on atoms. 



III. ENTANGLEMENT PURIFICATION OF ATOMIC STATES 



Entangled states are the most important resource for quantum information processing, including quantum tele- 
port at ion 



34 



37| . Generally, these applications 



35[, efficient quantum computation [36l and quantum cryptography 
work best with pure, maximally entangled quantum states. However, due to dissipation and decoherence, practically 
available states are likely to be non-maximally entangled, partially mixed, or both. To obtain useful entangled states, 
there are various entanglement purification schemes proposed, which differ in their purification range, the number of 
copies of the state they operate on. In this section, we focus on entanglement distillation which operate on a single 



copy, recurrence protocol which operates on two copies at each step. 



A. Bipartite Entanglement distillation 



Entanglement distillation 



i3, 
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12l | aims at preparing highly entangled states out of a supply of weakly entangled 
pairs, using local devices and classical communication only. The simplest protocol operate on a single copy of an 
arbitrary mixed state p and consist in the application of a POVM. Hence sequences of local operations are applied 
in such a way that for specific measurement outcomes the resulting state is more entangled than the initial state. 
Consider a special kind of mixed states including certain rank two states of the form 



/9 = i^|* + )(*+|+(l-F)|ll)(ll|, 

where = (|01)) + |10) /\/2 is one of the Bell states. We apply two POVMs (local flittering operators) 

A^B^\0) (0|+e|l)(l| 



(3.1) 



(3.2) 



on the mixed state and obtain the resulting state (non-normalized) 

p' = 1^-+) (*+| + (1 - F)e^ 111) (111 



(3.3) 
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The fidelity of the resulting state is given by 

F' ^Fe''/[Fe'' + {l~F)e\ (3.4) 

As expected, if e is small enough, the improved fidelity F' goes to 1 while the probability 

p = Fe^ + F)e^ (3.5) 

to obtain the Bell state 1^'"'') goes to zero. There is a tradeoff between the reachable fidelity of the output state and 
the success probability. 

Now we extend this method to any arbitrary bipartite mixed state and investigate how its entanglement changes 
under local operations and classical communications of the type p' ^ [A® B)p{A ® . Consider a pure entangled 
state |(/3) = a 1 00) + h\ll) (for simplicity, we assume (a > 6) G M) evolving through the same locally depolarizing 
channel Mi G {\/\ — pi, y/p/'ia^, ^J^pJ^tOy^ \fvT^'^z\ respectively, which is one of the most typical noisy channels in 
atomic systems. The resulting mixed state reads 



p = Y.'^Mi (g) Mj) \ip) {ip\ {M, (g) Mj)'^ 

id 
I 



(3.6) 



3a2(3-4p) + V 

6p - 4p2 

6p - 4p2 

\^ a6(3-4p)2 362(3 - 4p) +4p2 y 



a6(3 - 4p)2 





The optimal local filtering operations for a single copy of mixed state has been obtained by Verstraete et.al. in 
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12l |. Following their discussions, we introduce the real and linear parametrization of the p (R-picture) as 

p = Rij (Ji ® cr J, (3.7) 



% = Tiip{cri (g) aj)], 



(3.8) 



with (Jo — 1 and cri,cr2,o'3 arc the usual Pauli matrices. By using the singular value decomposition in Lorentz 
metric, the 4x4 matrix R can be decomposed as i? = LiSL^, where Li, L2 are some proper orthochronous Lorentz 
transformations, and S is of some normal form depending on R. The concurrence and the density matrix after the 



distillation are totally determined by E. In this case, we have E = diagjso, si, S2, S3} with 

so = ^(-V + 6p + A), 
2 

si = S2 = ga^(3 - 4p)2, 



S3 
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(V -6p + A), 



(3.9) 



The optimal local distillation operations can be easily obtained from Li and L2 directly and given by 

A®A*= tUiL^MT, (3.10) 



B®B*= THILIMT, 
^10 1 ^ 



(3.11) 



\/2 



V 



(3.12) 



110 

i -i 
10 0-1 

with M the usual Lorentz metric matrix M — diagjl, — 1, — 1, — 1}. After a straightforward algebra, the optimal 
POVM can be expressed as follows: 

/ \ / , 1/2 





A®B 



1 

1/2 



l-c 



1+c 







1+c 





(3.13) 



where 



12p + 9 + 2A 



(3.14) 



3(a2-62)(3-4p) 

B can be implemented with the setup in Fig. [T] by adding a phase shifter (QWP) behind partial polarizer, while A is 
obtained by a POVM in Eq. (2.10) followed by a single qubit rotation a^. on the atomic state by choosing the proper 

1 — c 



1+c 



The output state obtained after the application of local filtering operations can be represented as 

{A®B)p{A®B)'^ 



Tr 



{A®B)p{A®B)^ 



1 r,/ 



(3.15) 



with R' = E/sq. The distillation protocol produces the state p' with maximal possible entanglement of formation 



(EoF) according to the concurrence as introduced by Wootters 



23| 



r^t '\ rn -So + si + S2 - S3 
G(/9 ) = maxjO, 1 



max{0. 



2so 

a6(3 - 4p)2 + 4p2 - Qp 
-4p2 + 6p + A 



(3.16) 



compared to that of the initial state 



C{p) = max{0, ■ 



-So + Sl + S2 - 



= max{0, -[a6(3 - A.pf + V - 6p]}. 
9 



(3.17) 



Additionally the output state also has the maximum possible violation of the Clauser-Horne-Shimony-Holt (CHSH) 
version of inequality. Consider the expectation value /3p = Tr(/9S) for the operator 



(3.18) 



where (a, 6, c, d) are real unit vectors. In the CHSH inequality, the proposed value /?p, a combination of four polarization 
correlation probabilities, should not be more than 1 for local hidden variables theory (to obtain the value of the usual 
CHSH inequality, the factor 2 should be involved). In our case, the maximal CHSH violations for both initial and 
purified states are shown as 



Pp = ^a6(3-4p)2. 
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dp' = 



V2ab{3 - 4pj^ 
-4p2 + 6p + A' 



(3.19) 



In this paper we will focus on how to implement the entanglement distillation on neutral atoms, so next an- 



other typical noise channel — amplitude damping channel is considered. The corresponding super-operators are 

}. An entangled state \lp) = a\00) + b\ll) passes through the noise channel 



1 



M, G { 











y 



and then we obtain a mixed state 

+ b'^p'^ ab{l - p) 

b^{l-p)p 

b^{l-p)p 

y ab{l~p) 5^(1 -p)p J 

If the initial state is one of Bell state 1$"''). i.e. a = b = 1/V2, the 4x4 matrix R is shown as 



(3.20) 



R = 



^ I p ^ 

l-p 

p- 1 

p 2p2 - 2p + 1 y 



(3.21) 
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The optimal local filtering operation can be simplifies as 



A(g)B 



V 







1/2 X / 

» ' 



V 



V 



A 

1/2 





(3.22) 



Similarly, the POVMs A and B can be implemented by cavity-assistant interaction and linear optical elements with 

1/2 



the proper choices of e 



After the application of local filtering operations, the 4x4 matrix R' according to the output state p' is obtained 



as 



R' 



^f- + 1 - _p. 



1 









Vp^TT-p 

+ 1 _ p 





(3.23) 



with the concurrence 



C{p') = max{0, 



1-p 



(vV + i-p)' y 



-} 



(3.24) 



compared to that of the initial mixed state C(/9) = max{0, (1 — p)^}. The maximal CHSH violations for both initial 
and purified states are obtained as 

V2 



/3p = x/2(l-p),/3,. 



(3.25) 



Next we distill the entanglement from the general two-qubit mixed state. We let two qubits of the nonmaximally 
entangled states pass through the same phase damping channel in basis {|0) , |1)}, respectively. The corresponding 
super-operators are {-^/l — pi, y^(Tz}A ® {vT~^l, VP''^z}b- After decoherence, we obtain a mixed state shown as 

^ ' a6(-l-f 2p)2 ^ 







yafe(-l + 2p)2 








(3.26) 



For this kind of mixed states, the only local operation is an identity on the first particle and a POVM (unilateral 

h 



local filtering) on the second particle A _B = 1 < 



The POVM B can be implemented by cavity-assistant 



a 
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interaction and linear optical elements with the proper choices of e = a/b. After the application of local filtering 
operations, the output state is obtained as p' oc {A®B)p{A®By and in R-picture, the matrix R' can be decomposed 



as 



R' 



\ 



1 
1 
(l-2p)2 


The corresponding concurrences and CHSH violations can also be obtained and list as follows 

1 



(3.27) 



(1-2^)2^ 



C{p') 



2\ab\ 
1 

2\ab\ 



C{p) = (1 - 2pY 



/3p = v/1 + (1 - 2p)4. 



(3.28) 



The entanglement distillation protocol can be generalized to arbitrary two-qubit partially mixed state, since any 
nontrivial local operation and classical communication can be written as the form 7C/A U'j^^U-b U^, 

where J7a(b) and Wj^f^-^^ denote local unitary operations on the first (second) and 7 is a scale factor in the range 
< 7 < 1 and < a{(3) < 1. 

Figs. [5]and [3]show the concurrences and violations of CHSH inequality along with the parameter p for depolarizing 
channel and amplitude damping channel respectively. One can find that local filtering operations can always yield a 
new state with larger entanglement of formation. And there exist some cases where the initial mixed states do not 
violate the CHSH inequality while violate them after the local filtering operations. This also compatible with the 
results in 



101, 



11 



12l |. However, we can also find that although C{p) and /3p are both a entanglement measure, they 
are not just equal. Local filtering operations do not always lead to a larger violation of the CHSH inequality. In 
Fig. [21 if p is smaller than 0.16, the concurrence and violation of CHSH inequality of the state after the distillation 
process are greater than those of the mixed state before the distillation; if p is greater than 0.16, the concurrence of 
the state after the distillation is still greater than that of the initial mixed state, while the violation of the state after 
the distillation is even smaller. For the damping channel, it is similar and the critical point is p ~ 0.775. Thus the 
concurrence and violation of the CHSH inequality are two independent entanglement measures. 
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B. recurrence protocols 

In the following wc discuss another protocol to produce states arbitrarily close to a maximally entangled pure state 
by iterative application. That is so-called the recurrence protocol. 

Consider a mixed atomic state pi2 resulting from an imperfect distribution of |$^). We decompose p into two 
terms, p ~ p + pod, and express the density operator p in the Bell basis {|$^) , \^^) , |^^) , |^ )} and denote the 
diagonal elements in that basis by {A, B,C, D}. All off-diagonal elements in the Bell basis (pod) is made irrelevant 
by the protocol. Note that the first diagonal element A = ($+1 p |$^), which is the definition of the fidelity. 

Given two mixed states pi2 and pi'2', described by {A, B,C, D} and {A' , B' ,C' , D'} (and the irrelevant Pod) 
respectively, the following sequence of local operations obtains with a certain probability a state with higher fidelity 
and hence purifies the state: (i) application of a^^i (Ei Ox.i or 1 with probability 1/2 on p\2 and similar to pv2'\ 
(ii) the partial BSM on both p\2 and pv2' , i-C-, X]i=i.2 X]j=i,2 -^^/^ ^ /°i2Pi'2' , where M\ ~ and 
M2 = -P{|*+).[*-)}, which can be implemented with projection measurements shown in Sec II. We only keep the state 
P12 if * = J j i-c. the results in the final measiu'ement coincide in both states. 

The effect of (i) is to erase off-diagonal terms of the form which may contribute to the protocol. The 

operation in (ii) is to select states in atoms 1' and 2' which are eigenstates of Oz ® <Jx with eigenvalue +1 while 
eigenstates with eigenvalue -1 are discarded. In particular, we find that the remaining off-diagonal elements do 
not contribute and the action of the protocol can be described by the non-linear mapping of corresponding vector 
af — {A, B, C, D), ~x' — {A', B', C\ D'). The resulting state is of the form p + pod, and p will on average have diagonal 
elements given by 

A = (AA' + CC") /N, B = (BB' + DD') /N, (3.29) 
C = {BD' + DB') /N, D = {AC + CA') /N, 

and N = (A + D) {A' + D') + {B + C) {B' + C) is the probability of success of the protocol. This map is equivalent 



to the purification map obtained in Ref. [2j| for non-encoded Bell states. It follows that an iteration of the map- 
which corresponds to iteratively applying the purification procedure (i-ii) to two identical copies of states resulting 
from successful previous purification rounds-leads to a maximally entangled state. That is, the map has {1, 0, 0, 0} as 
attracting fixed point whenever A > B+C+D. We emphasize that all errors leading outside the subspace, independent 
of their probability of occurrence, can be corrected. The method is capable of purifying a collection of pairs in any 
state p, whose average fidelity with respect to at least one maximally entangled state with probability greater than 1/2. 
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The central part of recurrence protocol of entanglement purification — partial BSM can be implemented on neutral 
atoms via cavity-assistant interaction and linear optical elements shown above. 

Here we remark that the methods such as nested purification can also be applied via BSMs, which significantly 
reduced the required number nodes 25l |. 



IV. FEASIBILITY 

In this section, we analyze the feasibility of this scheme using numerical simulations in the experimental parameters. 

No particularly demanding assumptions have been made for the experimental parameters. The relevant cavity 
QED parameters for our system are assumed to be (507^,7) /27r = (27,2.4,2.6) MHz [29|. g'^/n'y ~ 117 ^ 1 places 
our system well into the strongly coupled regime. The cavity consists of two 1-mm-diam mirrors with 10 cm radii of 



curvature separated by 75 ^m |29j assuming the wavelength of the cavity mode is ~ 780 nm (the rubidium D2 line). 
The distance between two atoms d in an optical lattice is about 10 /im, which is larger than the waist of the cavity 
^ 5 /xm, leaving only one atom inside the cavity and its neighbor atoms outside for the gate operations. 

The evolution of the two-atom states is accomplished during the passing time of the single-photon pulse T ^ 
200/k — 10 ^s. The maximum velocity of the atoms in the transverse optical lattices and the maximum acceleration 
imparted are about 30 cm/s and 1.5g, respectively. Moving the proper atoms into and outside of the cavity is 
accomplished within a time tt ~ 100 /is. The gate preformation and the transport of atoms can be accomplished 
within the coherent time (dephasing) of the atoms, lasting milliseconds or hundreds of milliseconds, which depends 



on the sensitivity to the magnetic fluctuations of the internal atomic states 
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The partial polarizer has been implemented experimentally. For example, the partial polarizer with four slabs 



14| allows for a transmission probability Th = 0.18 (horizontal polarization) and Ty = 0.89 (vertical polarization) 
respectively. The transmission of horizontally polarized light is thus only 20% of that for the vertically polarized light, 
i.e. £ — 0.2. In principle, by tuning 0, the angle between the partial polarizer and the vertical axis, we can obtain 
any value of e (see Appendix for details). Hence, our scheme fits well the status of current experimental technology. 

Experimentally, major sources of deviations from the ideal setup are found to be: addressing errors, spontaneous 
emissions and long term interferometric phase instability. Since we only apply local operations, there is no interfer- 
ometer required in the setup of our proposal. Also, spontaneous emissions only lead to photon losses which can be 
detected, merely decrease the probability of success and thus have no contribution to a lack of fidelity. 

Some additional sources of fidelity degradation have been estimated by numerical simulations, which show that the 



14 



proposed scheme works remarkably well. In these simulations we considered the fidelity of the POVM shown in Fig. 
m which is the same to the fidelity of the output states p' obtained after entanglement distillation. Further decrease 
of fidelity is for example caused by shape mismatching between the input and output pulses. Fig. [2^ shows that we 
retain a high fidelity of the POVM under (^,7) /27r = (2.4,2.6) MHz (for all atoms with equal coupling coefficient 
g) and e = 0.2 —which means that J-povM is up to 99% for .g/27r > 27MHz [2^. Our method is also insensitive 
to randomness in the coupling rates caused by fiuctuations in the position of the atom, which can be numerically 
calculated by variation of g. As the figure shows, J^poym is extremely insensitive to this influence, as SJ-'pqvm, 
describing the change of the fidelity, stays below lO"'' for g varying from 27 MHz to 13.5 MHz. It is shown in Fig. [2] 
that the fidelity of POVM. 

From Sec. II (C), the fidelity of the operation of BSM is about Tbsm ~ -^povM- ^^^^ fidelity of the final state 
obtained after purification with the recurrence protocol is then AJ-^qm- 

The efficiency of this scheme is characterized through exact numerical simulations that incorporate various sources 
of noise and this demonstrates the practicality with the background of current experimental technology. 



In summary, we present a scheme to realize bipartite entanglement purification on neutral atoms via cavity-assistant 
interaction and linear optical elements. Two kinds of purification protocols are proposed here. That is the entangle- 
ment distillation which can be realized by POVM on atomic states and the recurrence protocol based on partial BSM. 
The proposal has been tested by numerical simulations including various sources of noise present in experiments, and 
was found to be feasible for experimental realization based on the status of current technology. 



In principle, we can always establish a relationship (shown in Fig. [5]) between the transmission probability and the 
angle 9 between the partial polarizer and the vertical axis based on Fresnel equation. At the upper interface, we 
obtain the transmission probabilities and as 



V. 



CONCLUSION 



Appendix 



2ricos6' 2cos0sin0' 



(Al.a) 



th = 



n'cosO + ncosO' sin(6' + 9') coa{e - 9') ' 
2ncos9 2cos6'sin0' 



71COS6' + n'cos9' sui{9 + 9') ' 



(Al.b) 
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where n and n' are the refractive indexes. Consider both interfaces, we obtain the total transmission probabilities as 



^ /IN sin20sin20' ,,„,n 
Tv = 4 = — 97 7- (A2.b) 



If substituting Eq. (A2.a) into (A2.b), we will obtain e = T\-^/T^ as functions of 6. 
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FIG. 1: (a) Schematic setup for the implementation of the POVM on atoms through photon-scattering. Moving the optical 
lattice takes single atoms in and out of the cavity. Scattered photon pulses leak out of the cavity after a reflection, pass throught 
the linear optical elements and are finally detected by two single-photon detectors. For a detection event in Dl, the outcome 
Pi is obtained directly; for a detection event in D2 we obtain Pi after a single-qubit rotation applied to the atom, (b) Relevant 
energy level structure of the atoms and their coupling configuration. 
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FIG. 2: The concurrence C (a) and violation of CHSH inequality /3 (b) along with the parameter of the depolarizing channel 
p. Here we show an example with the initial state 0.97 |00) + 0.25 The blue and red lines are for the mixed state before 
and after the distillation process. 




FIG. 3: The concurrence C (a) and violation of CHSH inequality j3 (b) along with the parameter of the damping channel p. 
Here we show an example with the initial state (|00) + |11)) /\/2- The blue and red lines are for the mixed state before and 
after the distillation process. 
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FIG. 4: (a) The fidelity of the POVM shown in Eq. (2.10) versus g/2-K with pulse duration T = lO^is, k/2-k = 2.4MHz, 
7/27r = 2.6MHz, e = 0.2. We have assumed a Gaussian shape for the input pulse with f{t) oc exp [— {t — T/2)^ / (T/S)^] . (b) 
It changes with g/ go- 




FIG. 5: The relationship between the transmission probabilities and the angle 9 can be always established based on Fresnel 
equation. 



